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We calculate the non-linear tunneling current through a molecule with two electron-accepting 
orbitals which interact with an intramolecular vibration. We investigate the interplay between 
Coulomb blockade and non-equilibrium vibration-assisted tunneling under the following assump- 
tions: (i) The Coulomb charging effect restricts the number of extra electrons to one. (ii) The 
orbitals are non-degenerate and couple asymmetrically to the vibration, (iii) The tunneling induces 
a non-equilibrium vibrational distribution; we compare with the opposite limit of strong relaxation 
of the vibration due to some dissipative environment. We find that a non-equilibrium feedback 
mechanism in the tunneling transitions generates strong negative differential conductance (NDC) in 
the model with two competing orbitals, whereas in a one-orbital model it leads only to weak NDC. 
In addition, we find another mechanism leading to weak NDC over a broader range of applied volt- 
ages. This pervasive effect is completely robust against strong relaxation of the vibrational energy. 
Importantly, the modulation of the electronic transport is based on an intramolecular asymmetry. 
We show that one can infer a non-equilibrium vibrational distribution when finding NDC under 
distinct gate- and bias-voltage conditions. In contrast, we demonstrate that any NDC effect in the 
one-orbital case is completely suppressed in the strong relaxation limit. 
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Optical spectroscopy of vibrational modes has pro- 
vided detailed information about the structure of 
molecules Q. In a similar fashion, tunneling-current 
spectroscopy @] may reveal microscopic details of single- 
molecule devices [E & IE IE HI • At low temperature 
electrons tunneling onto a molecular device excite dis- 
crete vibrations by spending some of their excess energy 
provided by the bias voltage. The resulting changes in 
the current with bias voltage are determined by the ef- 
fective potentials for the nuclear vibration in different 
charge states, in contrast to optical spectroscopy. This 
provides information on how the molecule is situated in 
a nanojunction (center of mass modes between the elec- 
trodes) Q and the role of different parts of the molecule 
in the transport (internal modes) [E IE ■ It is also of 
interest to control the electronic response of molecular 
devices by designing their mechanical properties through 
chemical synthesis. In this respect, the center of mass 
modes of a molecule in a nanojunction are less attrac- 
tive since they are difficult to tailor and most sensitive 
to relaxation through interaction with the surface of the 
junction electrodes Internal modes of a molecu- 

lar device, in contrast, are expected to remain more dis- 
crete. They reflect more intrinsic properties of the device 
amendable to chemical engineering. Recently, the cou- 
pling of such a tailored mode to the tunneling transport 
has been experimentally investigated in a C140 fullerene- 
dimer jfj. In view of the ongoing experimental efforts, 
it is of interest which current-spectroscopic features of a 
molecular device distinguish vibrational excitations from 
electronic ones. 

Theoretical works have already addressed several issues 
in the weak tunneling limit by focusing on a one-orbital 
model including the Coulomb blockade effect and cou- 
pling to a single localized vibration. At finite bias volt- 
age the electron tunneling will tend to drive the vibration 



out of equilibrium [U 13 El 13 Several mechanisms 
may reduce the accumulated vibrational ener gy: cou pling 
to the environment , the tunneling itself 0, Q, Hlf 
and also intramolecular vibrational energy redistribution 
due to anharmonic mode coupling. The effects of the 
renormalization of the tunneling by Franck-Condon (FC) 
factors on the vibrational distribution were discussed in 
detail [13|. These factors incorporate the effect of the 
change in the nuclear configuration when changing the 
electronic state and charge of the molecule. In the limit 
of strong electron- vibration coupling current suppression 
and a related super-poissonian current noise were pre- 
dicted as well as a weak NDC effect 0, M 

Strong 

NDC effects ^3 were found by assuming the electron- 
vibration coupling to increase with bias voltage. Finally, 
in the strong tunneling limit a Kondo effect due to charge 
(instead of spin) fluctuations was discussed for strong 
coupling to the vibration 0, 0] . 

In this paper we consider a molecule with two electron- 
accepting orbitals coupled to a single internal vibration 
with frequency lo. We consider the weak tunneling limit 
r <C T where T is the typical tunneling rate to the 
electrodes and T the electron temperature. The strong 
Coulomb charging effect is assumed to restrict the num- 
ber of electrons which can be added to the molecule to 
one. This introduces a competition between vibration- 
assisted tunneling processes associated with the different 
orbitals. We focus on the case where the orbitals are non- 
degenerate (splitting A) and asymmetrically coupled to 
the vibration (couplings Ai 7^ A2). Additionally, we as- 
sume that the vibrational distribution on the molecule 
can be driven out of equilibrium (no relaxation). For 
comparison we also discuss the technically simpler limit 
of strong relaxation. Due to the asymmetric coupling 
to the vibration the bare electronic splitting is renormal- 
ized to the observable splitting A which we will consider 
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as an effective parameter. We find that negative dif- 
ferential conductance (NDC) effects are enhanced com- 
pared with the one-orbital case and actually dominate the 
transport. Basically, the orbital coupled stronger to the 
vibration contributes little to the transport but functions 
as a trap due to a feedback mechanism in the tunneling 
transitions. The competing weakly coupled orbital gives 
the main contribution to the current but at the same 
time enhances the feedback by providing an additional 
path to the trap. Several types of features appear: (i) 
Franck-Condon progressions of alternating conductance 
resonances and anti-resonances resulting in current os- 
cillations which are robust against strong relaxation; (ii) 
anomalous current peaks of width oc T due to a strongly 
bias dependent redistribution of vibrational energy; (hi) 
isolated gate- and bias- voltage regions outside of which 
the current is strongly suppressed due to a stabilization 
of the charged state which couples stronger to the vibra- 
tion, (iv) when the excited orbital couples stronger to 
the vibration a voltage-controlled population inversion 
between the two charged electronic states takes place 
which is signalled by NDC occuring at special resonances. 
Importantly, the effects are based on the Coulomb block- 
ade and the intramolecular asymmetry (non-degenerate 
electronic states and asymmetric coupling to the vibra- 
tion). No asymmetric electronic wave function overlap 
with the electrodes needs to be assumed (similar to the 
NDC induced by to spontaneous emission [l9|l. This 
offers the interesting perspective of designing electronic 
transport properties of a molecular device by synthetic 
control of the electro-mechanical aspects. We note that 
for a related model of two degenerate orbitals current rec- 
tification was predicted assuming strong relaxation |2C| . 
There the asymmetric coupling to the electrodes is essen- 
tial. 

The paper is organized as follows. In Section [I] we in- 
troduce the two-orbital model and the master equations 
for the molecular state occupancies. We discuss the 
qualitative dependence of the FC-factors on the vibra- 
tional numbers and show that the Condon-parabola from 
optical spectroscopy is also a useful tool in transport- 
spectroscopy. In Section Hi Al we first discuss the NDC 
effect noted in Ref. fLj| for the one-orbital model and re- 
late it to the feedback mechanism. We then exhaustively 
discuss the effect of this mechanism for the two-orbital 
case in the limits where one orbital couples weakly to the 
vibration: A^ > 1 > \% ( Section ItTB|) and A^ < 1 < A2 
(Section 111 Cf) . In Section IIIDI we illustrate how the 
mechanism leads to more complex results for asymmetric 
strong coupling Af 5^ A| > 1. We present a comprehen- 
sive overview and discussion in Section IlIII 



I. MODEL 

Molecule. We consider the minimal model H = 
Hm + Tl r H r + Ht incorporating the molecule (M), 
the electrodes r = L, R and the tunneling (T) in units 
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The Hamiltonian Hm describes a molecule with two 
electron-accepting orbitals % = 1,2 (operators Ci ai en- 
ergies ei) with an energy splitting A = £2 — e±. Here 

rii = J2cr nicr > Ui ° ~ c ia c i<r- We assume throughout 
that double occupation of each orbital is completely sup- 
pressed due to strong local Coulomb interactions iij i.e. 
ni < 1. Simultaneous occupation of the two orbitals 
is similarly suppressed by the interaction v which intro- 
duces an important correlation: J^. iii < 1. Since Ui and 
v are the largest energy scales they do not enter explic- 
itly in any further way. The nuclear configuration of the 
molecule is assumed to be altered with respect to some 
coordinate Q when either orbital is occupied. The ef- 
fective potentials for the nuclear motion in these charged 
states (here in the harmonic approximation, frequency u) 
are shifted relative to the neutral state by i — 1, 2. 

The coordinate Q = (b + b^)/ \pl is normalized to the nu- 
clear zero-point motion by [Mlo)^ 1 / 2 (M = nuclear mass 
involved) and A; is dimcnsionless. Here excites the vi- 
brational mode by one quantum u> and P = (6— b j <)/\/2i 
is the conjugate nuclear momentum. The energy scale 
characterizing the electron-vibration coupling associated 
with orbital i = 1,2 is the change in the elastic energy 
at fixed nuclear configuration ujXf. Vibration-assisted 
processes are thus expected to lead to a progression of 
conductance resonances spread over a bias voltage range 
of at least ~ max{A|oj}. Typical energies uo of inter- 
nal vibrations observed experimentally range up to a few 
tens of meV 0, IS El • In general electron- vibration 
coupling is expected to be particularly strong for many- 
particle states of molecules which are characterized to 
a good approximation by occupation of a particular or- 
bital [2]]], which is typically the case for charged states 
of otherwise neutral molecules. Large relative displace- 
ments I Ai| > 1 of the nuclear potentials may be expected, 
for instance, when Q is an angle coordinate and the nu- 
clear configuration of the charged molecule is internally 
twisted relative to the neutral one (2^ . 
The electronic energy parameters in Eq. ^ are the rele- 
vant effective parameters for finite A^. These are related 
to their values for Ai — A2 — (indicated by a super- 
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— 2A1A2W. This is seen by diagonalizing the molec- 
ular Hamiltonian through a translation of the nuclear 



coordinates |23(, U = JJ, 
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The resulting 



Hamiltonian has the form £[J where \,i — 1,2 is elimi- 
nated. The electron becomes "dressed" with vibrational 
excitations (polaron) resulting in the renormalization of 
the energy parameters which we anticipated in writing 
Eq. ^ The renormalization of the charging energies is ir- 
relevant here since we assume them to be largest energy 
scales (i.e. uf ] > 2\X l \ 2 uj,i = 1,2 and > 2\\ 1 \ 2 \w). 
The correlations ri\ + ri2 < l,7i< < 1 are thus not weak- 
ened. Cases where a strong renormalization of the in- 
teraction becomes relevant were discussed in 0, Il8j . 
In contrast to the one-orbital model the renormaliza- 
tion of the orbital energies is important when the cou- 
pling to the vibration is asymmetric, Ai =/= A2. Then the 
electronic splitting is renormalized to an effective value 
A = £2 — ei = At ' + — A|) which can even have a 
different sign as A^ ' . Since only the excitation energy 
A is observable in the transport characteristics, we use 
it as an independent positive parameter i.e. the state 
1 by definition has the lowest renormalized energy e%. 
We are interested in the case where resonances related 
to orbital and vibrational excitations occur on the same 
voltage scale, i.e. ui ~ A < max{Afw}. The transport 
mechanism which we wish to illustrate operates in the 
limit of asymmetric coupling. This requires that either 
the lowest orbital couples strongly to the vibration or the 
excited orbital, see Fig.Q] We point out that in Hamilto- 
nian Q we have not written intramolecular terms which 
couple the two nuclear potentials 1 and 2. Such terms be- 
come important, for instance, when the electronic energy 
splitting is an integer multiple p of the vibrational ene rgy 
quantum, A = pu>. This has been discussed in Ref. |20j 
for the case A = and Ai = — A2. Here we avoid such 
degeneracies, i.e. we assume t <C min p {A — pa;}, where 
t is a tunneling amplitude between the electronic states. 
Furthermore, we can safely disregard electronic transi- 
tions induced by the nuclear motion near the crossing of 
the potentials 1 and 2 since for large asymmetric coupling 
the barrier separating the minima of potentials i = 1,2 
is (A/a>(Ai - A 2 ) -1 ± (Ai - A 2 )) 2 w/4 > uo. Below we will 
also present results for cases of moderate asymmetry of 
the vibrational couplings where such effects may start to 
play a role. These will serve as a simple starting point 
for the discussion of the strong asymmetry case and also 
illustrates the enhancement of NDC effects when multi- 
ple orbitals (instead of just one) are competing in the 
transport. 

The electrodes r = L,R are modeled by H r , Eq. J5J, 
as non-interacting quasi-particle reservoirs at electro- 
chemical potential fi r . The electrode-molecule tunneling 
Ht, Eq. 10, picks up the shift of the nuclear coordinate 
from the unitary transformation of the molecular oper- 
ators: H T = Efc«rr^4<xr e ~ A ' (b1 ^ b) ^ + h.c. Here c] CT 
creates a polaron state associated with the effective po- 
tential of electronic orbital i — 1,2. Since we consider 
here an intramolecular vibration we do not include a de- 
pendence of the bare tunneling matrix elements U r on 
the coordinate Q (shuttle-effect, cf. Ref. (HHU). 
Master equations. We are interested in the weak tun- 
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(a) A 2 «A! (b) A L «A 2 




FIG. 1: Schematic effective nuclear potentials for neutral (i = 
0) and charged (i — 1,2) electronic states of molecule, (a) 
Strongly coupled ground state (b) Strongly coupled excited 
state. 



neling regime, r <C T, where in addition the vibrational 
excitations can be resolved, T <C uj. We can describe 
the transport using diagonal density matrix elements 
P q (occupation probabilities, Yli=o Eg Pq = 1) j where 
q = 0, 1, 2, . . . is the vibrational number and i = 1, 2 de- 
notes the charged state with only orbital i occupied and 
the neutral state is labeled by i = 0. The transitions 
including the vibrational excitations q, q' are denoted by 
9 £3 l q i and 9 ^ 2 q i. The occupation probabilities are 
coupled by the stationary master equations 

i rq' 

Pi = =J2 (2WT^ 0? ,P° - W^^Pl) , (4) 

rq' 

where i = 1, 2, with transition rates (/ r (e) = (e^ e_Mr ^ T + 

ws q ^ q , = r« (l - / r (/4,_„)). (5) 

The addition energies for the transition i q > <— 9 are 

fi q ,_ q = e t + (q' -q)uj-aV g . (6) 

The gate voltage V g effectively varies [i relative to the 
ground-state transition energy /^J (a = capacitance ratio) 
and the bias voltage V > is applied symmetrically, 
Hl,r = /i ± V/2. The stationary current flowing out of 
reservoir r = L, R is given by {It + Ir = 0) 

!r = EE ( 2W i^0q'Pq> ~ W^P^) . (7) 
qq' i 

The equations for the one-orbital case with coupling Ai 
are obtained by simply discarding all P q , in Eqs. 
and (0 and are equivalent to those in Refs. 0, 0, 0, 
Il4l Il5j| . The current will change whenever a line in the 
(/1, V) plane is crossed corresponding to a right-electrode 
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resonance (j.r = fJ, q /_ q (positive slope in (/x, V) plane) or a 
left-electrode resonance hl — Mg'-« (negative slope). Im- 
portantly, due to the harmonic excitation spectrum only 
the change in vibrational number q' — q enters in the res- 
onance condition: transitions between all states i q i and 
q with fixed difference q' — q become allowed at a single 
resonance. A cascade of single-electron transitions can 
then lead to a significant population of high vibrational 
excitations, e.g. ig — > Oo — > i\ — > Oi — > 12 ~ * O2 — > ■ • • is 
a possible cascade for [xl > Mi > Mo > MR- 
Franck- Condon factors and Condon parabola. The tun- 
neling rates consist of two factors: T qq , — T n F q ,. The 
tunneling rates T™ = 27r|ii r | 2 p r between electrode r = 
L, i? (density of states p r ) and orbital i = 1, 2 determine 
the overall current scale. The FC factors i^i = F qql take 
into account that the stable nuclear geometry is changed 
when occupying orbital i: 

F qq ' = \(Q\XiW)\ 2 = e- X ^xf q - q ' 1 (L[^\X^) 2 ,(8) 

where L is the associated Laguerre-polynomial and q < 
q' . Note that the sign of A^ does not play a role. The gen- 
eral sum rule J2 q F qq , = ^2, q , F qq , = 1, guarantees that 
the current will saturate at large bias voltage to the value 
it would have without the vibrations ("electronic limit"). 
This holds only when the Xi do not depend strongly on 
the bias voltage (cf. 0). 

Without vibrations, asymmetry of the tunneling rates 
with respect to the orbital- and electrode- index gives rise 
to NDC l25| an d super-poissonian current noise [2(| , see 
also [2]], l28l l29| . Below we show that qualitatively differ- 
ent dependence of the FC factors F qq , on the vibrational 
numbers for state i = 1, 2 and the effective energy split- 
ting A give rise to NDC effects which can dominate the 
transport. We therefore set T lr = T,i = 1,2, r = L,R 
and restrict our discussion to V > since in this case 
I(—V) = —I(V). We symmetrize the stationary current 
/ = (It, — Ir)/2 and decompose it into a sum of pos- 
itive partial currents of the states weighted with their 
occupation: 

' = E« + ^EEw (9) 

q i q> 

7 " = r EE^'(/^' ? '- ? )-/«W- 9 )) (io) 

q' i 

4 = r E F «V (M/4'-*) - (11) 

For low T -C oj the partial currents are the FC factors 
F qq , summed over the transitions q i q > inside the 
bias window > A*g'- ? > M-R m Fig. [3 Note that the 
partial current of the neutral state I q has contributions 
from both charged states i = 1, 2 into which it can decay. 
One can understand the numerical results in almost all 
detail using the following simple graphical scheme. This 
approach works for multiple orbitals and also for a more 
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FIG. 2: Franck-Condon factors F* 1 in linear grey-scale 
(white = positive, gray=zero) for a transition between i q i 
and q involving a strong relative displacement Xi — 3. The 
Condon-parabola Eq. 1121 (white full line) separates the clas- 
sically forbidden and allowed regions. The partial current 
Il'/T, Eq. 1111 . is obtained by adding the FC-factors for tran- 
sitions in the bias window (strip of width and height V/ui be- 
tween the dashed white lines) horizontally for fixed q' . The 
partial current Iq/T, Eq. (11UI is obtained similarly by adding 
the FC-factors in this strip vertically for fixed q and adding 
the contributions from both states i = 1,2. 

general shape of the nuclear potential. Without going 
into the details, we comment on the basic points in the 
procedure. In Fig. [21 the FC factor associated with the 
transition between a pair of states q ^=5 i q > is depicted as 
function of q and q' . The change in the partial currents 
in (JSJ with increasing bias can be understood by drawing 
the bias voltage window in this figure. Only qualitative 
features of the FC factors are of importance which follow 
from simple quasi-classical arguments. The FC-factor 
F qq , is basically non-zero only in the classically allowed 
region delimited by the tilted Condon parabola 0: 



In this region the classical orbits of the nuclear motion in 
the shifted potentials intersect in phase-space. The FC- 
factor oscillates with q, q' with a quasi-classical envelope 
which varies algebraically on the scale of X\ [30|. The 
global maxima « 1 / A? are attained where the parabola 
touches the axes ( q ~ Xf,q' = and q — 0,q' » Af). 
In the classically forbidden region where the opposite of 
condition l|12|l holds the FC factor F qq , is exponentially 
small. In the forbidden regions where q' 3> q or q' <C q 
the nuclear momenta of the two motions are incompati- 
ble. These regions exist for any value of Ai . Starting from 
the allowed region the FC-factors eventually decrease ex- 
ponentially with increasing q' or q. On the other hand, 
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the forbidden region where q, q' <§C Xf is only well defined 
for A? ^> 1. In this case , initially increases expo- 
nentially with increasing q' or q <C A?. This "inverted" 
regime exists for nuclear potentials of more general shape 
with large relative displacements of the potential minima. 
As the bias window in Fig. [5] widens, the partial currents 
increase. For weak coupling Xf <C 1 the Condon-parabola 
is very narrow i.e. transitions which conserve the vibra- 
tional number have the largest amplitude. When the 
bias window reaches the vertex of this narrow parabola 
nearly all partial currents reach their maximal value at 
once. For strong coupling, Af 3> 1, the parabola is very 
broad and the partial currents show a slow exponential 
increase as the bias window widens. The complex trans- 
port characteristics of the two orbital model considered 
here follow from two intramolecular asymmetries between 
the orbitals: (1) the bias window covers different parts 
of the Condon-parabola due to the electronic splitting A 
and (2) the Condon-parabolas are qualitatively different 
due to asymmetric coupling to the vibration. 
Strong relaxation. We now prove an important restriction 
on the occurrence of NDC in the limit where the vibra- 
tional excitations completely relax before each tunneling 
event due to some dissipative environment. This limit im- 
plies the factorization ansatz P % q — P % P q with the vibra- 
tional equilibrium distribution P q = e~ qu l T (\ — e~"/ T ). 
We can then reduce the equations 10} to an effective elec- 
tronic three-level problem with effective bias-voltage de- 
pendent rates (cf. [Tol [T^ l obtained by averaging over 
the equilibrium distribution: 

WSLi = £ HVv p </ > Wfc-o = E W l^o g > p q' (13) 

qq' qq' 

(these vary monotonically with V) and Wo<— i = 
Y,r W o^i >W»<-o = Y,r W i^o- Tne stationary current 
in electrode r = L,R reads ( II + Ir = 0): 

£,:Ko-^) 

In Appendix^ we show that in this case for two orbitals 
the current can be reduced by increasing the bias voltage 
at resonances related to the left electrode, fiL = = 
1, . . . (for V > 0), i.e. lines with positive slope in the 
(/i, V) plane. Here the transition i q +k <— g becomes al- 
lowed. Any NDC along a resonance line with positive 
slope is thus a proof of a non- equilibrium vibrational dis- 
tribution on the molecule. For one orbital the current 
can never decrease with V in this limit. Finally, we con- 
sider the limit where in addition to the strong relaxation, 
the transitions g ^ i q >, i = 1,2 are not correlated i.e. 
the renormalized Coulomb interaction [TtI is zero: 
v' = 0. It is readily shown ( Appendix lB|l that in this case 
the current increases monotonically with bias voltage for 
any number of orbitals. The strong Coulomb correlations 
are thus essential for NDC effects. 

Intermediate relaxation. Basically all the physics is cap- 
tured by considering the opposite limits of negligible and 



strong relaxation. We have confirmed this by consider- 
ing intermediate regimes where we add a relaxation term 
J2 q i Wq^q' P % q , to the right-hand side of the equation for 

PL Eq. J3J. We considered an environment 11] with 
either ohmic (s = 1) or sub-ohmic (s = 0) spectral func- 
tion J(E) = -f\E/io\ s for which the relaxation rates are 
Wq^q. = J(u>(q - q'))[±N{u(q - q'))} for q ^ q' where 
N(E) = {eP E — 1) _1 . We briefly discuss the results for in- 
termediate relaxation when interesting deviations from a 
simple interpolation between non-equilibrium and strong 
relaxation limit occur. 



II. RESULTS 

We now present results for the stationary current / 
(Eq. 0) and differential conductance dl/dV for symmet- 
ric tunneling rates r ir = T,r = L,R,i = 1,2. Through- 
out the paper we set the temperature to T = 0.025w. 
Gray-scale plots of dI/dV(fi, V) have been given different 
linear scale factors for dl/dV ^ to clarify the voltage 
conditions under which NDC occurs. The NDC mag- 
nitude can be inferred from the presented /(V) curves. 
We note that the results for a molecule with two neutral 
states and one charged state are simply obtained by in- 
verting the sign of /i — and modifying the discussion 
accordingly. 

A. Strongly coupled single orbital - Feedback 
mechanism and weak NDC 

In Ref. ^3 the one-orbital model in the limit of strong 
coupling to the vibration (Af ^> 1) was shown to exhibit 
a current suppression which was related to "avalanches" 
of vibration assisted tunneling processes which also leads 
to super-poissonian current noise effects 0, 0] . Addi- 
tionally, for asymmetric gate voltages a weak NDC effect 
was noted in the absence of relaxation. Here we focus on 
this weak NDC effect and additional small current peaks 
(not discussed in Ref. ^{|) which are visible in Fig. |3| 
as black lines and white-black double lines, respectively. 
The mechanism responsible for this is based on an en- 
ergy asymmetry (induced by the gate potential) and the 
qualitative features of the FC-factors. This mechanism 
will also play a role in the transport involving two (or 
more) orbitals where it results in much stronger effects. 
We therefore consider this simple case in some detail. 

Current suppression. The differential conductance plot- 
ted in Fig. |3 is symmetric about /i = /xj up to unim- 
portant differences in amplitude due to spin degeneracy 
of the orbital (without which we would have exact sym- 
metry) so we only discuss the case \i > /ij; the opposite 
case follows from interchanging the roles of the electronic 
states <-» 1. For low voltages V <C 2ajXf states with 
low vibrational number are predominantly occupied and 
the current is exponentially sup pre ssed both in the limit 
of weak and strong relaxation {l5l ITr3 |. This is related 
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FIG. 3: One-orbital model with strong electron-vibration 
coupling Ai = 5.0. Differential conductance in gray-scale 
(gray: dl/dV = 0, white / black: dl/dV ^ 0) as function 
of bias V and gate voltage jj,. 
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FIG. 4: Current-voltage characteristics for increasing fi — 
fil = Otv, 0.25iv, 1.75oj (going down) in Fig [3] for non- 
equilibrium vibrations (Eq. J7J). Inset: small current peak 
for fi — fig — 0.25a; which for larger /i — /iq develops into a 
current drop. 



to the classically forbidden region q, q' <C \\ in Fig. [21 
where the FC factors depend exponentially on the vi- 
brational numbers. In the limit of strong relaxation the 
current suppression is simply due to the exponentially 
small partial currents of the few vibrational states which 
are thermally occupied at low T <C lo. In the absence of 
relaxation the FC-factors also prevent the excited states 
from actually becoming occupied for V <C 2uiXf . The re- 
sulting non-equilibrium vibrational distribution induced 
by the tunneling is "equilibrium-like" as was noted be- 
fore 0, . This is due a type of feedback mechanism 
in the tunneling transitions, as we will now explain. In 
principle at finite bias voltage arbitrarily high vibrational 
excitations can be accessed via cascades of single-electron 




FIG. 5: Relevant transitions for Fig. For clarity only 
vibrational energy differences are indicated. For fi > fil a 
feedback cascade of transitions keeps the state lo nearly fully 
occupied leading to NDC and current peaks. 



tunneling processes. However, at low bias voltage tran- 
sitions which lie outside the bias window in Fig. [21 (i.e. 
q — > l q ' resp. q <— l q i ) correspond to large changes 
of the vibrational energy and have exponentially larger 
rates. Once the initial states for these transitions start 
to be occupied, the total rate for populating the low- 
est states (transitions outside the bias window) becomes 
much larger than the total rates of its decay (transitions 
inside bias window) . As a result only the low- lying states 
are occupied due to the large asymmetry between the 
rates. Compared with the strong relaxation limit vibra- 
tional excitations are slightly more favored and therefore 
the current suppression is less severe [l5| in this limit 
(not shown). The central observation is that although 
the occupations decrease strongly with vibrational num- 
ber this is compensated by the exponential increase of the 
partial currents. The main contribution to the current 
comes from the excited states. The "inverted" depen- 
dence of the FC-factors on energy (vibrational numbers) 
thus stabilizes the lowest vibrational state and enhances 
the sensitivity of the current to the small occupations 
of the vibrational excitations. This is at the basis of the 
weak NDC and small current peaks which we will discuss 
now. 

Weak NDC and current peak. Apart from the quantita- 
tive effect on the current suppression, there appear inter- 
esting qualitative new features in the absence of relax- 
ation. The many visible resonance lines in Fig. Inform a 
Franck- Condon progression which extends beyond volt- 
ages ~ 2u>Xi. Interestingly, for asymmetric gate energy 
/i — (Aq > u>/2 the current can be reduced at the res- 
onances hl — /ij.,fc = 1,2,... (black features on right 
in Fig. EI). Here the transitions q — > l q +k 7 q = 0, 1, . . . 
become allowed and one could naively expect the cur- 
rent to increase since this favors the population of the 
excited states l q+ k,q > which are responsible for the 
main current contributions. Actually, the opposite hap- 
pens. Due to the asymmetric gate-energy [i — /^J > a 
feedback mechanism is active which involves a cascade of 
single-electron transitions. This is illustrated in Fig. [51 
In the regime of voltages where NDC occurs the charged 
vibrational ground state lo is stabilized, not only relative 
to the neutral states (due to the Coulomb blockade) but 
also relative to the vibrationally excited charged states 
l q ,q > 0. Due to the asymmetric gate energy state lo 
has less transitions which depopulate it than which pop- 
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ulate it, see Fig|5| Due to the strong increase of the 
FC-factors with vibrational number this asymmetry in 
the rates causes a nearly complete occupation of lo for 
asymmetric gate energies. Upon increasing the bias, at 
resonances where the excitations l q +k become accessible 
other excitations g ' with q' w q + k are also favored via 
subsequent tunneling processes involving small changes 
in the vibrational number. These subsequent transitions 
are allowed already at low V. The excitations g > decay 
with large rates back to the the state lo and its first few 
excitations, since they change the vibrational number by 
a large amount. Importantly, the reverse of the latter 
transition, g / «— lo is not allowed at low T, see Fig. |S] 
Therefore the occupation of lo is effectively increased at 
the expense of the excited states which contribute most 
to the current and NDC occurs. 

Small current peaks of width oc T occur in the interme- 
diate region < /i — /ij < w/2 (white-black double lines 
in Fig. and inset of Fig. . These signal a redistri- 
bution of the vibrational energy when the bias is tuned 
through the resonance. In this case the above feedback 
mechanism can only become effective when the transition 
energy lies sufficiently close to hl . Thus initially the cur- 
rent rises but once the excited states become sufficiently 
populated they start to relax via the feedback and the 
current drops again. 

Finally, Fig. 0] shows that even though the absolute 
current-step amplitude exponentially increases with in- 
creasing bias voltage 0], the NDC becomes relatively 
less pronounced. Careful inspection reveals that at suf- 
ficiently large gate energy |/i — the resonances ini- 
tially correspond to current drops but with increasing 
voltage turn into peaks and finally become current steps. 
The increasing bias eventually compensates for the gate- 
asymmetry and the feedback mechanism becomes inef- 
fective. 

In the limit of strong relaxation at low T -C uj only tran- 
sitions inside the bias window along the q — and q = 
axis in Fig. [3 play a role (cf. Eg. 1 1 3|> - The charged vi- 
brational ground state is stabilized due to a gate voltage 
(i — /LtJ > and the strong relaxation. However, no NDC 
or current peaks can occur at the resonances discussed 
above since the vibrational distribution is not affected by 
the tunneling in this limit. The feedback mechanism is 
cut off: after each single-electron tunneling process the 
excitation relaxes on a much shorter time scale and the 
next tunneling process starts from a vibrational ground- 
state again. Indeed one can show explicitly that in the 
limit of complete relaxation the NDC and the current 
peaks disappear in the one-orbital model for arbitrary 
spin and orbital degeneracy (see Appendix IX)). We note 
that at resonance lines [ir = fj,_ k for k = 1,2,... the 
current always increases, independent of the relaxation. 
Here the transitions l q — > q +k become allowed whereby 
lo can decay and repopulate the excited states which 
carry the current. 
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FIG. 6: Differential conductance dI/dV{n,V) for Ai = 1.4, 
A2 = 0.1, A/lo — 2.5. Current steps (white lines with nega- 
tive slope) due to stronger coupled orbital 1 become current 
drops (dark lines) as soon as the weakly coupled orbital 2 
contributes to the current. 

B. Strongly coupled ground state 

We now demonstrate that in a two-orbital model with 
asymmetric coupling to the vibration already at moder- 
ate coupling a new but weak NDC effect occurs which 
is robust against strong relaxation, in contrast to the 
weak NDC in the one-orbital model. Additionally, the 
non-equilibrium feedback mechanism responsible for the 
weak NDC in the one-orbital model leads to strong NDC 
effects for two competing orbitals. These general state- 
ments carry over to the case of multiple non-degenerate 
orbitals with asymmetric coupling occupied by at most 
one electron due to Coulomb blockade. We start our dis- 
cussion with an intermediate case where only the weak 
NDC occurs. 

1. \{ > 1 ^> A! Weak NDC - Current oscillations 

The differential conductance and typical I-V curve in 
Fig. and [7| respectively, display a number of features 
which can be understood by considering the two orbitals 
invidually. Vibrational excitations of the stronger cou- 
pled state 1 form a FC progression of resonances far 
beyond V ~ 2X\uj. In contrast, the excitations of the 
weakly coupled orbital 2 only show up as a single dl/dV 
peak at \xj, = The resonances /il = [i\ of the 

stronger coupled state 1, correspond to current steps for 
k ^ [A/u>]. Interestingly, once the weakly coupled state 
2 has started to contribute to the current, \xl > Moj these 
turn into anti-resonances, k > [A/uj], (dark lines in Fig-IHl 
with negative slope). The resonances, /j,r = n_u,k > 0, 
always correspond to current steps (white lines in Fig. 
with positive slope). A distinctive feature is that the 
current steps and drops due to the excitations of the 
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FIG. 7: Current-voltage characteristic for (j, — (i,q = 0.375lj 
in Fig. HJ for non-equilibrium (full, Eq. Q) and equilibrium 
vibrations (dashed, Eq. 1141 1. The oscillatory pattern of cur- 
rent steps and drops sets in only when orbital 2 contributes 
to the current (here V > 5.0cu). 



stronger coupled state have opposite gate voltage depen- 
dence. Since they are of the same order of magnitude 
they give rise to current oscillations on the slowly satu- 
rating background. We note that for identical parameters 
the one-orbital model (i.e. with either Ai =0.1 or 1.4) 
produces negligible NDC effects. The origin of the en- 
hancement of NDC is the following. For moderate gate 
energy /i — [1q the partial currents of states l q are small 
compared with those of states 2 q and q , cf. Fig. |21 At 
resonance lines /il = Mfc the transitions l q +k g be- 
come allowed, and the occupations of the states with the 
larger partial currents are reduced due to the Coulomb 
correlation between the two orbitals. The asymmetry 
between the partial currents is only present when the 
weakly coupled orbital 2 is accessible, therefore the cur- 
rent drops at these resonances only once both orbitals 
are accessible, i.e. /ig ^ /zr- In contrast, the res- 

onances fiR = fJ,_j,,k > correspond to current steps 
since here states l q are depopulated in favor of the states 
with larger partial currents. 

Although the current oscillations in Fig.[B]seem very sim- 
ilar to those for the one-orbital model in Fig. [21 there is 
an important difference: here the NDC is not completely 
suppressed in the strong relaxation limit, although larger 
values of the dominant coupling Ai are required for vis- 
ibility comparable with the limit of no relaxation. In 
Appendix^ we prove that in the strong relaxation limit 
a drop of the current can only occur along resonance lines 
Mi = Ml (negative slope) if it occurs. A condition for the 
visibility of NDC is that the total current below the res- 
onance is larger than compared with the partial current 
of the orbital causing it. This requires a weakly coupled 
orbital, A2 -C 1, with large partial currents in combina- 
tion with a stronger coupled orbital, Af > 1. These are 
roughly the same conditions as for the visibility of the os- 
cillating current in the non-equilibrium limit. However, 



we point out that under identical bias and gate voltage 
conditions the NDC need not to be visible in both the 
non-equilibrium and equilibrium limit, see for instance 
Fig. Hand Fig. HUI below. 

In summary, the current oscillation occurs due to the 
competition between transport channels with signifi- 
cantly different partial currents. It is a Coulomb repul- 
sion effect: the current calculated without the effective 
correlations (v = 0), always increases with V (see Ap- 
pendix [HJ. 



2. Af > 1 > Al Strong NDC 

When the charged ground state couples strongly to the 
vibration , Af 3> 1, and the electronic excitation lies low, 
A < u>, we find the typical structure of Fig. EI The cur- 
rent oscillations are clearly visible again and this set of 
resonances needs no further discussion. An obvious dif- 
ference with Fig. is the finite gap > A for any gate 
voltage. In fact, in addition to the Coulomb blockade 
regime where /iL < Mo or Mo < Mfl resp., the current is 
suppressed in the entire strip /j,q > Mz, > Men M-R > Mo (i- e - 
where the excited orbital 2 is not yet accessible). This is 
due to the exponentially small FC-factors , ,q,q'<ti Af 
of the lowest orbital which now couples strongly to the 
vibration (cf. Sec. Ill Ajl . When increasing the bias volt- 
age above the gap the resulting current depends strongly 
on the order in which additional transitions become al- 
lowed i.e. on the gate voltage. Basically four different 
situations can occur which are labeled (a)-(d) in Fig. 
and the relevant transitions are schematically indicated 
in Fig. El 

(a) Stabilization of the charged ground state. In the 
upper-right region in Fig. [HI the current may be expected 
to flow: the transitions 2 q > => g are energetically al- 
lowed for at least q = q' = 0, for which the FC-factor 
is large, Fq 1. Instead the current is strongly sup- 
pressed and increases in small steps of increasing height 
with increasing bias. This is very similar to the situation 
in Fig. 0] where the transport is dominated by a single 
orbital with Af ^> 1. Between states l q ' and q the feed- 
back mechanism discussed in Sec. Ill Al is operative which 
keeps state lo almost fully occupied with increasing bias, 
see Fig. I^a) and compare with Fig. |SJ The presence of 
orbital 2 further enhances the feedback since the states 
q (which feed back into lo) can now also be populated 
via cascades of transitions involving orbital 2. 

(b,c) Isolated region. The small diamond-shaped region 
in Fig. |S1 at finite voltage A < V < 2uj is remarkable. In- 
side it the current is non-zero and beyond any of its four 
defining boundaries in the plane of gate- and bias volt- 
age the current is completely suppressed. A typical I{V) 
curve through this region is shown in Fig. ^3 The fact 
that the region is bounded by a NDC line with positive 
slope proves that it must be caused by non-equilibrium vi- 
brational effects (SecQJl, see also the strong relaxation re- 
sult in Fig. 1101 This region can be reproduced by truncat- 
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FIG. 8: dI/dV{n,V) for Ai = 5, A 2 = 0.1 and A/w = 0.5. 
The isolate region at low bias is due to the small excitation 
energy A < uj and disappears for A > uj. The current is 
strongly suppressed once the weakly coupled state 2 can be 
excited, /il > fli and hr < fi 2 _ 1 (thick black lines). A cur- 
rent peak (white-black double line) appears when this orbital 
becomes accessible at /il = Ho- 



(d) (a) 




(c) (b) 




FIG. 9: Relevant transitions in low bias section of Fig. |H1 
Indicate are transitions which cause NDC to occur when in- 
creasing the bias along (a)-(d). For clarity only vibrational 
energies are indicated. Arrow thicknesses indicate relative 
magnitude of the transition rates between vibrational excita- 
tions of a pair of states q ^ i q i, q, q' = 0, 1, 2, . . . where i = 1 
or i = 2. Rates between different pairs of states are of differ- 
ent order of magnitude since different couplings (Ai or A2) is 
involved. The cascades of transitions give rise to a feedback 
into the lowest vibrational state lo of strongly of the coupled 
orbital. 



ing the spectrum to 5 states: g , 2 q , q = 0, 1 and lo, see 
also Fig. [5J At the low bias side, Coulomb blockade and 
the small FC factor Fq discussed above are responsible 
for the current suppression. Inside this region only the 
transitions 2o 0o ^ lo are allowed. The stationary oc- 
cupations follow from Eq. ®: P ° = 1/5, Pq 1 = P 2 = 2/5 
and the current is: 

1 « + « = §r (15) 
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FIG. 10: /(V) for fi-fil = 0.375lj in Fig.[H][case (b)], for non- 
equilibrium (full, Eq. J7J) and equilibrium vibrations (dashed, 
Eq. 1141 1. The isolated current plateau at low bias V < 2u> 
and the subsequent current suppression are non-equilibrium 
effects. For this gate energy both the current oscillations and 
the rising background current are enhanced when one allows 
for strong relaxation (in contrast to the one-orbital model). 



This is less than 2r/3 and 4r/5, the maximal current 
through one and two orbitals (without the vibration), re- 
spectively, due to the partial occupation of the strongly 
coupled state lo with suppressed partial currents. Note 
that state lo is not yet blocking the transport. At the 
high bias side, the current becomes suppressed when the 
first neutral excited state 0i can be reached either via 
the cascade 0o — > 2i — > 0i (NDC line with negative slope, 
case (b) in Figs. Elands or 2 — > 0i (NDC line with pos- 
itive slope, case (c)). State 0i can decay to lo when an 
electron enters the molecule through either junction i.e. 
the reverse transition 0i <— lo is suppressed at low tem- 
perature T <C uj — A. Now state 1 is almost fully occu- 
pied since it is populated much faster than it can decay: 
it is blocking the transport since the current it limited 
by the very small sum of its decay rates. The feedback 
loop thus involves both the weakly and strongly coupled 
state: the competition between the two orbitals which 
couple asymmetrically to the vibration causes the NDC 
to be much stronger compared with the one-orbital case 
(Sec lII A} . From the above it follows that the diamond- 
shaped region disappears for excitation energy A > lo. 
However, the feature (a) and (d), which we discuss now, 
will still be present. 

(d) Current peak. At the resonance line /i^ = /Xq where 
the weakly coupled orbital 2 starts to participate a sin- 
gle large current step could be expected. Remarkably, 
NDC occurs in the middle of this resonance, produc- 
ing a current peak (white-black double line with negative 
slope in Fig. |SJ) whose width is proportional to the elec- 
tron temperature T. This sharp features stands out be- 
tween the thermally broadened plateaus in Fig. ^2 The 
origin of the peak is a strong competition between the 
charged states 1q and 2q in the narrow energy window 
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FIG. 11: I(V) for n - /4 = -1.125w in Fig. |U [case (d)] 
for non-equilibrium (full, Eq. Q) and equilibrium vibrations 
(dashed, Eq. 1141 '). The current sets on with a peak of width 
oc T, the other local maxima are thermally smeared plateaus. 
The downward trend of the equilibrium curve only applies at 
low voltages V < 2Aioj, it saturates to the "electronic limit" 
at larger voltages. 



\fj,L — (J,q\ *~ T, involving the feedback via neutral excited 
states q . This is most simply illustrated by considering 
gate energies A/2 — uj < \i — \i\ < A/2 — uj/2 where 
a minimal set of 6 states is sufficient to understand the 
peak, see Fig. Eld). At the rising side of the peak the 
rate of the transition 2o «— 0o (through the V depen- 
dence in the Fermi-function, cf. Eq. |SJ has increased 
sufficiently to enhance the current relative to the very 
small value supported only by the transitions to / from 
the blocking state, 1q Oo- Due to the gate energy 
the simplest feedback loop involves a cascade of 6 tran- 
sitions: Oo — > 2 — > 0i — > 2i — > O2 — > lo- This feedback 
initially remains ineffective since the excited state O2 is 
not yet sufficiently populated. The vibrational distribu- 
tion is equilibrium-like and the current follows the result 
for equilibrated vibrations, see Fig. ^2 As one increases 
V the occupations of the excited states in the feedback 
loop increases and a maximum is reached. Here the feed- 
back dynamically starts to trap the molecule in the state 
lo- The occupations of the excited states and the current 
now start to decrease and reach a lower value (although 
higher than before the peak). The current peak thus 
signals this redistribution of the vibrational energy in a 
small bias window. We note that the current is not com- 
pletely suppressed: this only happens beyond the second 
strong NDC line /il — Mi (°f- Fig- Hlf > as may be un- 
derstood by considering Fig.^Jd). The peak can thus be 
considered a precursor of the full onset of the feedback 
mechanism. For lower < A/2-w a similar argu- 

ment involving more than 6 states explains why the peak 
becomes a step and why simultaneously the strong NDC 
along — fif is further enhanced. 

Intermediate relaxation. Upon increasing the vibrational 
relaxation rate 7 (cf. SecP) starting from zero, the strong 



relaxation result is approached as expected. The relax- 
ation cuts off the cascade of transitions leading to the 
blocking state and reduces the importance of the feed- 
back for the transport. Now the NDC becomes more 
pronounced at resonances fii, = fj,^, k = 1, 2, ... where the 
transitions q — > l q +k become allowed. These enhance 
the occupation of lo due to relaxation and suppress the 
current. However, this approach is rather slow at low 
bias voltages: the NDC lines marking the isolated region 
remain clearly visible. 

In summary, the strong NDC lines = [i\ and /j,r = 
y?_y are associated with excitations of the state with 
weakly coupled to the vibration (state 2). However, the 
strongly coupled state 1 is actually blocking the trans- 
port. The weakly coupled state allows an excess vibra- 
tional energy to accumulate on the molecule (through 
a cascade of tunneling processes) which is subsequently 
spent to trap the molecule in the strongly coupled state 
(in a single tunneling process). The blocking state can 
thus be reached under very general energetic conditions. 
Therefore NDC effects become strong when 2 (or more) 
orbitals which couple asymmetrically to the vibration 
compete in the transport. Finally we note that Fig. [S] 
is reminiscent of the signatures of spin-blockade of tun- 
neling. There the resonance line marking the transition 
between ground states can be terminated at finite bias 
and the current is only recovered when excited states be- 
come accessible [3lJ. The resonance line thus shows a 
kink. Here the kink in the resonance line is more drastic 
since also transitions to many vibrationally excited states 
are suppressed. Such details are of importance to distin- 
guish NDC due to spin excitations in molecules from 
effects due to vibrational excitations. 



C. Stronger coupled excited state 

We now consider the case opposite to Section III Bl 
where the excited orbital 2 takes up the role of the block- 
ing orbital due to a strong(er) coupling to the vibration 
and the lower orbital 1 is weakly coupled. Vibration- 
assisted tunneling processes now stabilize the electroni- 
cally excited state of the charged molecule i.e. an excess 
charge and energy may be stored on the molecule by the 
feedback mechanism. We focus on the resulting qualita- 
tive differences with respect to Sec. Ill HI resulting from 
this population-inversion controlled by the bias voltage. 
As a simple reference point, we start again from an in- 
termediate case. 



1. A? <C 1 < A2 Weak NDC - current oscillations 

Fig. E| shows the typical differential conductance for 
A > uj. Similar to Fig. the weakly coupled state 
(state 1) basically shows up as one large current step and 
the many resonance lines correspond to the excitations 
of the strongly coupled state (here state 2). However, at 
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FIG. 12: dI/dV(n, V) for Ai = 0.1, A 2 = 1.4, and A/lu = 2.5. 
The oscillations are similar to those in Fig ll2l The current 
drops (black lines with negative slope) are again due to the 
stronger coupled orbital 2. 

the resonances [ir, = \j? k , k > [A /lu] we have current drops 
for small |/x — /x§ | (i.e. away from the charge degeneracy 
point) and current steps for large \fi — [Aq\. Therefore, 
around the charge degeneracy point fi = /ij the first few 
excitations beyond the electronic excited line [ij, = [i\ 
show up as current steps and only at higher bias voltage 
turn into current drops, in contrast to the case where the 
lower orbital couples stronger to the vibration (Fig. EJ. 
Up to now we have only found features in the differen- 
tial conductance (either positive or negative) when in- 
coming electrons excite the vibration with their excess 
energy. This is expected at low temperatures T <C lu. 
Interestingly, in Fig. ^| a small current step at the reso- 
nance fiL — M-i with negative slope is visible which cor- 
responds to absorption of the vibrational energy by an 
incoming electron despite the low temperature T <C lu 
and moderate bias. The slightly enhanced current may 
be understood as an effect of significant heating of the 
molecule by the vibration assisted-tunneling. 

2. A? < K A! Strong NDC 

When the coupling to the charged excited state be- 
comes strong, A| 3> 1 the case of most interest is that 
of a higher excited orbital, lu < A <C \\uu. The typi- 
cal structure is shown in Fig. In contrast to Fig. [S] 
we do not have a gap here since the change in the nu- 
clear configuration of the two ground states is now small: 
the current starts to flow at the edges of the Coulomb 
blockade region [it, = /ij and /j,q = [ir. Two NDC lines 
stand out in Fig. 1131 (dark) where the current is signifi- 
cantly suppressed: /jl = /i\ and hr = Again, the 
appearance of the latter NDC line with positive slope 
is a proof that non-equilibrium vibrations play a role 
(Sect. 0. At these lines the transitions l g +x <— q and 
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FIG. 13: dI/dV(n,V) for Ai = 0.1, A 2 = 5.0, and A/w = 
2.5. Similar to Fig. |H| a two strong NDC lines define a low 
bias structure and a current peak, which now occurs at fi > 
fig . The progression of NDC lines extending up to large bias 
voltage is again due to the strongly coupled excited orbital 2. 
Remarkably, this NDC already sets in at low V, before this 
orbital can be directly accessed: this signals absorption of 
vibrational energy by incoming electrons. 



(d) (a) 




FIG. 14: Same diagram as Fig.|5|but now for low-bias section 
of Fig. 1131 The feedback mechanism now involves more exci- 
tations Q q ,l q ,q < A/uj > 1. Note the state which is stabilized 
by this non-equilibrium mechanism , 2o, is an electronically 
excited state i.e. we have a population inversion. 



lq — > 0q+i between the neutral and the weakly coupled 
charged state become allowed which enhance the occu- 
pation of the excited states Q q . The latter feed back to 
the strongly coupled state 2q with rates which increase 
exponentially with q and the electronic excited orbital is 
predominantly occupied. We thus have a bias-controlled 
population inversion between ground- and excited state 
of the charged molecule due to their asymmetric coupling 
to the vibration. We discuss the four different situations 
labeled (a)-(d) in Fig. ^| for which the relevant transi- 
tions are schematically indicated in Fig. ^] 
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(a) Current peak. A current peak appears along the line 
M-R = Mo f° r M — Mo > w /2 (white-black double line on 
the right in Fig. I18[l . Compared with Fig. El this peak 
occurs on the opposite side and remains visible up to 
much higher values of ^ — /xj. The mechanism causing 
the current peak is analogous to that discussed in Section 
III B 21 with the roles of 1 and 2 interchanged. However, 
more vibrational excitations are involved in the feedback 
mechanism depicted in Fig. O since A > ui. For moder- 
ate gate energy \i — the current peak is a precursor to 
the population inversion: the feedback is only completely 
activated beyond the resonance above it, ^jlr = 
There the transitions l q — > q +i become allowed and 
the strongly coupled state 2o is predominantly occupied 
suppressing the current (upper-right region in Fig. I13|l . 
For large asymmetric gate energy /x — /xj 3> oj/2 the peak 
actually marks the population inversion: the NDC at the 
peak gains in amplitude at the expense of the NDC above 
it at [i\ = (opposite to Fig. (SJ. In this case state 
2o can already be reached by many feedback cascades 
Oo —* lfc — * 2 , 1 < k < [(jJt— (Xq)/u)], once the escape from 
orbital 1 (Oo <— lo) becomes possible. Therefore the pop- 
ulation inversion and current suppression are complete at 
the peak. For a low- lying excited orbital, A < to, this is 
in fact the general situation since the cascade of transi- 
tions involved in the feedback is shorter. 
(b,c) Isolated region. For \fi — < uj/2 and V < 2lu 
we have an isolated region in the sense that the current 
reaches a local maximum value (diamond shaped region 
at bottom of Fig. 1 1 H|> . This region does not occur for 
A < uj (opposite to Sec. Ill B2l where A > oj suppresses 
the isolated region). Within this region only the transi- 
tions 1q ^ are allowed, Pq = 1/3, Pq = 2/3, and the 
current 

r « Jjitf + i/iitf = §r (16) 

equals the maximum current which a single orbital (with- 
out the vibration) can carry. When going along (b) and 
(c) in Fig. 1131 we next cross the resonance lines (11 = Mi 
and /ir = respectively and the current decreases 
as discussed above. The current suppression is complete 
once both transitions g — > l q +i and l q — » 9 +i are al- 
lowed. 

(d) Absorption by incoming electrons. In Fig. 1131 res- 
onances /x^ = M_ = 1>2 (negative slope) are vis- 
ible which correspond to absorption of the vibrational 
energy by an incoming electron (2 q <— q+ k). What is 
remarkable compared to Fig. 1121 is that the current de- 
creases here. This is another signature of the popula- 
tion inversion due to the feedback mechanism. Vibra- 
tional energy can accumulate on the molecule through 
previous sequences of tunneling events involving only the 
weakly coupled orbital 1. The molecule is then "brought 
to a standstill" when, in a single tunneling process, an 
electron with an energy deficit matching the total ac- 
cumulated vibrational energy enters: 2o <— Ofc. For 
/x| < /i£ < /ig there are [A/a;] such resonance lines where 
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FIG. 15: I(V) for /x-xx£ = 0.0 in Fig^ror non-equilibrium 
(full, Eq. 0) and equilibrium vibrations (dashed, Eq. (114B ') . 
The current plateau is caused by the weakly coupled ground 
state 1. The current suppression is due to the strongly cou- 
pled excited state 2 at A = 2.5a; which is only accessible by 
a single-electron tunneling process for V > 5.0ui (cf. Fig. |HJ. 
However, a cascade of such processes, whereby vibrational en- 
ergy accumulates on the molecule, allows this state to become 
occupied already at low V, Figs. I14f b1 .Cel. 

such a new trapping process becomes possible. Note that 
this can not be understood as an effect of heating since 
at these resonances the occupation of vibrational excited 
states is suppressed due to the feedback and the current 
is reduced. 

Intermediate relaxation. Compared with Section III B 21 
the non-equilibrium effects are more sensitive to relax- 
ation since here roughly 2 A/a; tunneling events comprise 
the feedback mechanism instead of 2 (cf. Figs. Inland 1 14| . 
This restricts the minimal vibrational quality factor Q 
(writing the vibrational relaxation rate as 7 = w/Q, cf. 
Sec. [J for the observation of effects due the feedback 
mechanism to Q > 2 A/a; for Af <S 1 <S A|. In con- 
trast, for 3> 1 3> A| the requirement is Q > 2. For 
the cases discussed here only Q ^> 5 is required. This 
is confirmed by our calculations for intermediate values 
of the relaxation rate 7. Interestingly, when increas- 
ing the vibrational relaxation rate 7 starting from zero, 
the dependence of the amplitude NDC line with positive 
slope [1r = p\ (proof of a non-equilibrium distribution) is 
non-monotonic: it is initially weakened and then regains 
amplitude and remains clearly visible with increasing 7. 
Also the current peak shifts to larger values of /1 — /ij and 
V but remains visible. The low-bias effects of the excited 
state are however suppressed since cascades responsible 
for the population inversion effect are cut off by the re- 
laxation. 

In summary, it is remarkable that in all discussed cases 
(a)-(d) the excited state 2o dominates the current at 
low bias where it can not be reached from the neutral 
ground state 0o by a single-electron tunneling process. 
The strong deviation from equilibrium is induced by cas- 
cades of single electron tunneling processes. 
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FIG. 16: dI/dV{p, V) for Ai = 5.0, A 2 = 1.1 and A/lu = 0.5. 
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-3-2-10123 
FIG. 17: dI/dV(fi, V) for Ai = 5.0, A 2 = 1.1 and A/lu = 2.5. 

D. Non-degenerate strongly coupled states 

Having analyzed the above cases in detail, we can now 
restrict ourselves to a brief classification of the results for 
non-degenerate states A > where both Af ^ A| > 1. 
Here the feedback mechanism discussed above produces 
more complex results. The basic change is that when the 
weakest coupling is increased the feedback mechanism be- 
comes less efficient at populating vibrational excitations 
of the neutral molecule (compare with the feedback mech- 
anism for a single orbital Sec. Ill A|) . More vibrational 
excitations of the weakly coupled orbital and the neutral 
state must be accessible by a single tunneling process in 
order to fully activate the feedback and trap the molecule 
in strongest coupled orbital. The patterns of NDC lines 
will therefore extend over a broader range of applied volt- 
ages. Indeed, a glance at Figs. 1161191 already shows 
that more NDC lines are visible. Also, there are more 
NDC lines with positive slope, which proves that the de- 
viations from an equilibrium vibrational distribution are 
stronger. For A^ ^> A2 > 1 the NDC effects are strongest 
for the case A < u> presented in Fig. EH Compared with 
Fig. [S] the isolated region defined by strong NDC lines 
is repeated a number of times to the left and it extends 
further to the right. Also, the current peak at the reso- 
nance line fiL = has shifted further to the left. The 
extended /(V) plateaus have a width fixed by u> — A, 
independent of the gate voltage (compare with the weak 
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FIG. 18: dI/dV{pL, V) for Ai = 1.1, A 2 = 5.0 and A/lu = 0.5. 
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FIG. 19: dI/dV(fi, V) for Ai = 1.1, A 2 = 5.0 and A/lo = 2.5. 

NDC in Sec. Ill B II where the current steps and drops 
have opposite gate- voltage dependence). For the case 
A > lj presented in Fig. El the low bias structure disap- 
pears and the current peak along fix, = Mo becomes the 
dominant feature. For A2 3> > 1 the NDC effects are 
strongest for the case A > lj presented in Fig. The 
two strong NDC lines in Fig. ^| have developed into a 
"checkerboard" pattern of such lines. These correspond 
to excitations of the weaker coupled state. In addition 
more resonances due to the strongly coupled state ap- 
pear. For A < 10 the current is more suppressed at posi- 
tive gate energies, Fig. ^| 

III. SUMMARY AND DISCUSSION 

We have calculated the non-linear current through a 
molecule with two non-degenerate electron-accepting or- 
bitals coupled asymmetrically to an internal vibration in 
the limit of weak tunneling to the electrodes. We found 
that due to the interplay of Coulomb blockade and non- 
equilibrium vibration-assisted tunneling NDC effects be- 
come amplified and pervasive in comparison with a one- 
orbital model. The only resonances where we consistently 
find current steps correspond to an electron tunneling 
off the molecule starting from the charged state coupled 
strongly to the vibration. At all other resonance lines the 
dl/dV may become negative depending on the electron- 
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vibration couplings and applied voltages. A weak and 
strong NDC effect may be distinguished, which require 
the larger of the two electron-vibration couplings to be 
moderate and strong respectively. 

The weak NDC effect is found at resonance lines where 
an electron can tunnel onto the molecule resulting in the 
charged state coupled stronger to the vibration. This 
effect only occurs when two (or more) orbitals are com- 
peting in the transport. The current steps and drops 
occur at bias positions with an opposite gate-voltage de- 
pendence and give rise to current oscillations. We proved 
that this type of NDC is robust against strong relaxation 
of the vibrational distribution on the molecule due to a 
dissipative environment. Any NDC at other resonances 
conditions is a proof of a non-equilibrium vibrational dis- 
tribution on the molecule. The current oscillations may 
even become amplified by strong relaxation depending on 
the applied voltages. 

Strong NDC effects appear at the first few resonance lines 
associated with the state weakly coupled to the vibra- 
tion. This is a non-equilibrium effect which is typically 
weakened by relaxation processes. Cascades of single- 
electron tunneling processes involving the vibrational ex- 
citations of the weakly coupled state provide a feedback 
which rapidly populates the strongly coupled state. The 
latter thus acts as a blocking state which is almost fully 
occupied. The few, strong NDC lines correspond to the 
activation of the feedback mechanism and can have the 
same gate voltage dependence as the current steps, in 
contrast to the weak NDC above. In a one-orbital model 
the feedback mechanism is also active but produces only a 
weak effect due to the absence of a competing orbital. An 
anomalous current peak of width oc T appears when the 
feedback mechanism becomes effective only sufficiently 
deep inside a resonance. The peak signals the crossover 
of the vibrational distribution from equilibrium to non- 
equilibrium. 

Interestingly, the blocking state can be the vibrational 
ground state of either charged state whichever couples 
stronger to the vibration. When the electronic ex- 
cited orbital couples most strongly the NDC signals 
a voltage-controlled population inversion between the 
charged states induced by the vibration-assisted tunnel- 
ing. Also, new resonances appear associated with an elec- 
tron entering the molecule and absorbing vibrational en- 
ergy stored on it (despite the low temperature) where the 
current is suppressed. 

The NDC effects are due to asymmetry of the orbital 
energies and couplings to the internal vibration which 
are intrinsic properties of the molecule. One can thus 
tailor the electronic response of the device by molecu- 
lar engineering. In contrast to other NDC effects, we 
do not require detailed assumptions of orbital- and/or 
electrode- specific electronic wave-function overlap with 
the electrodes 0, |25[ |23> l34l nor bias- voltage dependent 
coupling to the vibration |12j . 

For the interpretation of transport spectroscopy exper- 
iments on molecular devices an important result of our 



work is that multiple orbitals may be relevant for effects 
at voltages where only a single orbital would seem to mat- 
ter. For instance, an excited orbital may already domi- 
nate the transport by cascades of single tunneling events 
at a low voltage where this state is not directly accessi- 
ble from the neutral ground state, cf. Fig. 1131 Similarly 
the charged ground state may completely dominate the 
transport even at a high bias voltage where a far "better 
conducting" excited orbital is already directly accessible 
cf. Fig. 

We have used a basic parameterization of the nuclear 
potential surface of the electronic states and some com- 
ments are appropriate. For one, the nuclear potential 
shape may be anharmonic in the coordinate Q consid- 
ered here. The resulting qualitative changes may be de- 
termined from the FC-factors for these potentials, when 
plotted in similar fashion as Fig. The main results 
are not sensitive to the fine details of these factors but 
only to their large-scale dependence on the vibrational 
numbers due to the shift of the potentials, which can be 
established by quasi-classical considerations. Secondly, 
anharmonic terms in the nuclear potential may also cou- 
ple the mode Q to other internal modes which we have 
not considered here. When this coupling is strong for a 
large number of such other modes or a Fermi- (nonlinear- 
resonance is involved, intramolecular vibrational energy 
redistribution may relax the vibration. Generally, this 
will become more important for large molecules. The ef- 
fect of relaxation has been discussed. If however, only 
a few other modes couple strongly to Q, say one, an in- 
teresting two mode problem occurs. A treatment of the 
effects of such multi-mode dynamics pll | on the tunneling 
transport, lies outside the scope of the present paper. 
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APPENDIX A: STRONG RELAXATION 

We consider the limit r <C 7 <C T where the vibra- 
tional excitations have completely relaxed before each 
tunneling event due to a dissipative environment. The 
tunneling rates T lr are not assumed to be symmet- 
ric. With the factorization ansatz P % q — P l P q , P q = 
e -qu/TQ _ e ~u/T^ we can rec i uce equations (@J to an 
effective electronic three-level problem with voltage de- 
pendent rates (|13f) . The stationary probabilities and the 
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current can be explicitly given (r = R,L): 



g; =2 g^o ;P o = 



1 + 



(Al) 
(A2) 



We can now find a simple explicit condition for the oc- 
currence of current steps or drops with increasing bias 
voltage. Consider an increase of the positive bias V — > V 
such that one additional transition involving orbital i 
comes into the bias window through a resonance with 
electrode r = L, R, \i r = for some k — 0, ±1,±2, . . . 
(for V < interchange L «-> i? below). For simplicity we 
consider the values of the current on the two subsequent 
plateaus: only two transition rates are then changed, 
Wi<_o W/<-o an d Wo<-i — * Wo<-i the changes being 
related by <5W Vi = -e-^tfWT^ (cf. Eqs.EHSJ. The 
change in the stationary current 5I r = I' r — I r may be 
calculated at either electrode r = L,R from 114|) (since 



I'r + r, 



It 



0): 



Ir Wo^ Wo^i) \ I r 



(A3) 



Here f — R,L denotes the electrode opposite tor — L, R. 
In order to have NDC at a resonance where /il becomes 
larger than [i\ we require 51 l /II < 0, which, using r = L 
in Eq. |13J|, gives W£_i - h < 0. The rate of es- 
cape through junction R at voltage V (including the 
new transition in increased bias window V') must thus 
be smaller than the current II at initial voltage V . It 
is readily seen that this condition cannot be fulfilled in 
the case of only one orbital: Wq < _ 1 — I L = Wr^_ 1 +Ir = 
2Wft P° + WqLiC 1 ~ pl ) > °- However, for two (or 
more) orbitals it is possible to satisfy this requirement. 
To see if NDC may occur at a resonance where fj,R drops 
below [A we use r = R in Eq (1 A3j) . leading to the require- 
ment Wq^—Ir = Wq^+Il < 0. which can not be satis- 
fied for any applied voltages since It, > for V > 0. The 
current must increase at such resonances. Thus for fully 
equilibrated vibrations NDC can only occur at resonances 
related to the left electrode, = fi — fx\,k = 0,..., 



where electrons can enter the molecule by an additional 
tunneling process i q+ k <— q - These resonances corre- 
spond to lines with positive slope in the (fi, V) plane. 
Any NDC occurring at an other resonance is a proof of a 
non-equilibrium vibrational distribution. This proof can 
be trivially extended to iV orbitals correlated by Coulomb 
charging (maximally one extra electron). It also does not 
depend on the FC-factors involved, although the ampli- 
tude of the possible NDC may be small for a particular 
choice. 

APPENDIX B: UNCORRELATED 
VIBRATION-ASSISTED TUNNELING 



We consider the special case where the renormalization 
of the interaction (due to the polaron effect) compensates 
the Coulomb repulsion effects, i.e. v = — 2ljAiA 2 = 0. 
Now we have to include the doubly charged state of the 
molecule n\ — l,ni — 1 (di-anion). We denote the di- 
agonal density matrix elements by P™ 1 " 2 , where q is the 
vibrational number and nin2 denotes an electronic state 
with occupations n, = 0,1 of orbital i — 1,2. The occu- 
pations are coupled by the stationary master equations 

(cf. H3) 



E 



^w 0q ^i q + (i- ni )J2 2 w iq 



pn 1 n 2 
1 



+ J2[^iW lq ^ 0q 'P q r 2 + 2n 2 W 2q ^ 0q >P q/ l0 
+ (1 - n^Wo^P^ 2 + (1 - n 2 )Wo q ^ 2q 'P q , 1 ' 



nxn.2 



1. In the case where 



together with J2 ni n 2 P 
the vibration is assumed to be completely equilibrated, 
pnin 2 _ p™i"2 the kinetic equations can be decoupled 
into equations for the occupations of two uncorrelated 
"channels", P l = £ ni „ 2 S lnt P nin2 with the rates (15): 
P* = = 2Wi^ (l - P l ) ~ Wo^iP 1 . The current is then 
a sum of independent contributions of the individual or- 



bitals: I r = J2i In where I r i — (Wq+ 



2WL 



For a single orbital the current monotonically increases 
in the limit of strong relaxation (Appendix ^) and the 
same thus holds for two (or more) uncorrelated orbitals. 



[1] G. Herzberg, Molecular spectra and molecular structure 

(Van Nostrand, 1950). 
[2] L. Kouwenhoven, C. Marcus, P. McEuen, S. Tarucha, 

R. Westervelt, and N. Wingreen, Mesoscopic Electron 

Transport (Kluwer, 1997), p. 105. 
[3] H. Park, J. Park, A. K. L. Lim, E. H. Anderson, A. P. 

Alivisatos, and P. L. McEuen, Nature 407, 52 (2000). 
[4] J. Park, A. N. Pasupathy, J. I. Goldsmith, C. Chang, 

Y. Yaish, J. R. Petta, M. Rinkoski, J. P. Sethna, H. D. 

Abruna, P. L. McEuen, et al., Nature 417, 722 (2002). 
[5] A. N. Pasupathy, J. Park, C. Chang, A. V. Soldatov, 

S. Lebedkin, R. C. Bialczak, J. E. Grose, L. A. K. Donev, 



J. P. Sethna, D. C. Ralph, et al., Nano Lett. 5 (2005). 
[6] L. H. Yu, Z. K. Keane, J. W. Ciszek, L. Cheng, M. P. 

Stewart, J. M. Tour, and D. Natelson, Phys. Rev. Lett. 

93, 266802 (2004). 
[7] L. H. Yu and D. Natelson, Nanotechnology 15, S517 

(2004). 

[8] L. H. Yu and D. Natelson, Nano Lett. 4, 79 (2004). 

[9] S. Braig and K. Flensberg, Phys. Rev. B 68, 205324 

(2003) . 

[10] S. Braig and K. Flensberg, Phys. Rev. B 70, 085317 

(2004) . 

[11] D. Boese and H. Schoeller, Eur. Phys. Lett. 54, 668 



(2001). 

[12] K. D. McCarthy, N. Prokofev, and M. T. Tuominen, 

Phys. Rev. B 67, 245415 (2003). 
[13] A. Mitra, I. Aleiner, and A. J. Millis, Phys. Rev. B 69, 

245302 (2004). 

[14] J. Koch, M. Semmelhack, F. von Oppen, and A. Nitzan, 

cond-mat/0504095. 
[15] J. Koch and F. von Oppen, Phys. Rev. Lett. 94, 206804 

(2005). 

[16] J. Koch, M. Raikh, and F. von Oppen, cond- 
mat/0501065. 

[17] P. S. Cornaglia, H. Ness, and D. R. Grempel, PRL 93, 
147201 (2004). 

[18] P. S. Cornaglia, H. Ness, and D. R. Grempel, cond- 
mat/0409021. 

[19] M. H. Hettler, W. Wenzel, M. R. Wegewijs, and 
H. Schoeller, Phys. Rev. Lett. 90, 076805 (2003). 

[20] G. A. Kaat and K. Flensberg, cond-mat/0411173. 

[21] H. Koeppel, W. Domcke, and L. S. Cederbaum, Adv. 
Chem. Phys. 57, 59 (1984). 

[22] M. Cizek, M. Thoss, and W. Domcke, cond-mat/0411064. 

[23] I. G. Lang and Y. A. Firsov, Sov. Phys. JETP 16, 1301 
(1963). 

[24] J. Koch, F. von Oppen, Y. Oreg, and E. Sela, cond- 



16 



mat/0405453. 

[25] M. Hettler, H. Schoeller, and W. Wenzel, Nano-Physics 
and Bio-Electronics: A New Odyssey (Elsevier, Amster- 
dam, 2002), pp. 303-322. 

[26] A. Thielmann, M. H. Hettler, J. K6nig, and G. Schon, 
cond-mat/0406647 (2004). 

[27] A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. Lett. 
92, 206801 (2004). 

[28] A. Cottet, W. Belzig, and C. Bruder, Phys. Rev. B 70, 
115315 (2004). 

[29] W. Belzig, Phys. Rev. B 71, 161301(R) (2005). 

[30] A. Perelomov, Generalized Coherent States and Their 
Applications (Springer- Verlag, Berlin, 1986). 

[31] D. Weinmann, W. Hausler, and B. Kramer, Phys. Rev. 
Lett. 74, 984 (1995). 

[32] C. Romeike, M. R. Wegewijs, and H. Schoeller, cond- 
mat/0502091. 

[33] M. Hettler, H. Schoeller, and W. Wenzel, Eur. Phys. Lett. 
57, 571 (2002). 

[34] A. Thielmann, M. H. Hettler, J. K6nig, and G. Schon, 

Phys. Rev. B 68, 115105 (2003). 
[35] M. R. Wegewijs and Y. V. Nazarov, Phys. Rev. B 60, 

14318 (1999). 



